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Practical applications of chaos require the chaotic orbit to be robust, defined by the absence of periodic 
windows and coexisting attractors in some neighborhood of the parameter space. We show that robust 
chaos can occur in piecewise smooth systems and obtain the conditions of its occurrence. We illustrate this 
phenomenon with a practical example from electrical engineering. 



It has been proposed to make practical use of chaos in 
(S| communication [Q, in enhancing mixing in chemical pro- 
cesses Q and in spreading the spectrum of switch-mode 
power suppies to avoid electromagnetic interference ||, [| . 
^ ^ In such applications it will be necessary to obtain reliable 
Operation in the chaotic mode. 

. It is known that for most smooth chaotic systems (take 
CO the logistic map [9 for example), there is a dense set of pe- 
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r iodic windows for any range of parameter values. There- 



on fore in practical systems working in chaotic mode, slight 
inadvertent fluctuation of a parameter may take the sys- 
H tern out of chaos. The question is, how to guarantee that 
there is no periodic window for a given range of param- 
i eter values and the maximal Lyapunov exponent remains 
^ positive throughout the range? In this Letter, we address 
^ this problem. 

O ■ We say a chaotic attr actor is robust if, for its param- 
os* eter values there exists a neighborhood in the parameter 
• space with no periodic attractor and the chaotic attrac- 
rN tor is unique in that neighborhood. It is known that ro- 
^| bust chaos cannot occur in smooth systems. In this Let- 
ter we show that such situations can occur in piecewise 
smooth maps and obtain the conditions of existence of ro- 
bust chaos. 

We first give a practical example from electrical engi- 
neering to demonstrate robust chaos. The circuit shown in 
Fig.|l| is known as the boost converter. It consists of a con- 
trolled switch 5, an uncontrolled switch D, an inductor L, 
a capacitor C and a load resistor R. When the controlled 
switch is turned on, the current in the inductor increases 
and energy is stored in it. When the controlled switch is 
turned off, the stored energy in the inductor drops and the 
polarity of the inductor voltage changes so that it adds to 
the input voltage. The voltage across the inductor and 
the input voltage together "boosts" the output voltage to 
a value higher than the input voltage. Such circuits are 
widely used in regulated dc switch- mode power supplies. 

Regulation of the output current is achieved by con- 
trolling the switching by current feedback — known as 
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Figure 1: The current mode controlled boost converter 



"current-mode control". In this control logic, the switch 
is turned on by clock pulses that are spaced T seconds 
apart. When the switch is closed, the inductor current 
increases till it reaches the specified reference value I re f- 
The switch opens when i = I re f . Any clock pulse arriv- 
ing during the on period is ignored. Once the switch has 
opened, the next clock pulse causes it to close. 

We obtain a discrete-time model by observing the state 
variables at every clock instant. There are two ways in 
which a state can evolve from one clock instant to the next. 
If the on-time T on = L(I re f — i n )/Vi n is less than T, the 
evolution between observation instants includes one on pe- 
riod and one off period. Since the clock period typically is 
much smaller than the characteristic time of the LCR cir- 
cuit, we assume the waveforms to be linear between clock 
instants. By neglecting the higher order Taylor terms, the 
two dimensional map for T on <T is derived as: 
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On the other hand, if the clock pulse arrives while 
i < Iref, the switch remains on between the observation 
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instants. If T on > T, then the map takes the form 
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The borderline between the two cases is given by the 
case where the current reaches I re f exactly at the arrival 
of the next clock pulse, i.e., I border = I re f — V in TjL. The 
resulting map, therefore, is piecewise smooth. 
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Figure 2: The bifurcation diagram and the Lyapunov spec- 
trum of the boost converter. The parameter values are: 
c = 220/iF, I ref = 0.5A, V in = 30V, T = 400/is, L = 0.1H 

The bifurcation diagram and the Lyapunov spectrum of 
the boost converter is presented in Fig.||. It may be noted 
that there is no periodic window or coexisting attractor in 
the parameter range R= [241, 500] f2. The chaotic attrac- 
tor therefore satisfies the conditions of robustness. 

We now obtain the general conditions of occurrence of 
robust chaos. Let f(x,y;p) be a two-dimensional piece- 
wise smooth map which depends on a single parameter p. 
Let T p , given by x = h(y, p) denote a smooth curve that 
divides the phase plane into two regions R a and R^. The 
map is given by 
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x,y e R a , 
x,y e R b 
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It is assumed that the functions fa and fa are both con- 
tinuous and have continuous derivatives. The map / is 
continuous but its derivative is discontinuous at the line 
T p , called the "border". It is further assumed that the 
one-sided partial derivatives at the border are finite. We 



study the bifurcations of this system as the parameter p 
is varied. 

If a bifurcation occurs when the fixed point of the map is 
in one of the smooth regions R a or R^, it will be one of the 
"standard" types, namely period doubling, saddle-node or 
Hopf bifurcation. But if the bifurcation occurs when the 
fixed point is on the border, there is a discontinuous change 
in the elements of the Jacobian matrix as p is varied. A 
rich variety of bifurcations have been reported (|, |?], ||] 
in this situation, which have been called border collision 
bifurcation. We show that under certain conditions border 
collision bifurcation results in robust chaos. 

It has been shown |(J that by a change of coordinates, 
any piecewise smooth map can be reduced to the normal 
form (H) in some small neighborhood of the fixed point 
undergoing border collision bifurcation. 



Gn 



TL 




/ X 


Sl 


! 




tr 




(' 


Sr 




V y 




, for x < 0, 
, for x > 0, 
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where x and y are the new coordinates for which the border 
is along the line x = 0, dividing the phase space into two 
halves L and R. p is the new parameter, which is obtained 
by scaling p. tl and Sl are the trace and determinant of 
the Jacobian matrix in side L. tr and Sr are the corre- 
sponding values in side R. Since the trace and determinant 
are invariant under change of coordinates, these quantities 
are the same as in map /, calculated in the neighborhood 
of the point where border collision occurs. As the parame- 
ter p is varied through zero, local bifurcations depend only 
on the values of tl^Sl^tr, and Sr appearing in (||) and 
therefore it suffices to study the bifurcations in the normal 
form (||) in exploring the border collision bifurcations in 
the piecewise smooth map (|l|). 

The fixed points of the system in the two sides are given 

by 
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and their stability is determined by the eigenvalues \\^ = 
l ' 







48). 



It may be noted that if 

t l >(1 + S l ) and t r <(1 + Sr) 



(3) 



then there is no fixed point for p < and there are two fixed 
points, one each in L and i?, for p > 0. The two fixed points 
are born on the border at p = 0. We call this a border 
collision pair bifurcation. An analogous situation occurs 
if tl < (1 + Sl) and tr > (1 + Sr) as p is reduced through 
zero. Due to symmetry of the two cases, we consider only 
the parameter region (H). 
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If (1 + Sr) > r R > -(1 + Sr) then for /i > 0, the fixed 
point in L is a regular saddle and the one in R is an attrac- 
tor. This is like a saddle- node bifurcation occurring on the 
border. Since this region in the parameter space always 
has a periodic attr actor for /i > 0, we exclude this region 
from our analysis when looking for chaotic behavior. The 
condition tr = — (1 + 5r), results in a nongeneric situation 
where all points on the line joining the points ( i+s R •> 0) 

and (0, — Y+§^ ) are fixed points of the second iterate. We 
therefore concentrate on the parameter space region 

t l >(1 + Sl) and TR <-(1 + 5 R ) (4) 

and investigate the property of the attractor for fi>0. We 
first consider the case 1 > Sl > and 1 > Sr > 0. 

For (||), L* is a regular saddle and R* is a flip saddle. 
Let Ul and Sl be the unstable and stable manifold of L* 
and Ur and Sr be the unstable and stable manifold of 
R* respectively. For (||), all intersections of the unstable 
manifolds with x = map to the line y = 0. Since one 
linear map changes to another linear map across the x = 
line, Ul and Ur experience folds along the x-axis. And 
all images of fold points will be fold points. By a similar 
argument we conclude that Sl and Sr fold along the y- 
axis, and all pre-images of fold points are fold points. 

Let Ail, ^2L be the eigenvalues in side L and Ai#, 
X2R be the eigenvalues at side R. For condition (||), 
Ail > A2L > and > Xir > X2R. The stable eigenvec- 
tor at R* has a slope mi = (—Sr/Xir) and the unstable 
eigenvector has a slope = (—Sr/X2r). Since points 
on an eigenvector map to points on the same eigenvector 
and since points on the y-axis map to the x-axis, we con- 
clude that points of Ur to the left of y-axis map to points 
above x-axis. From this we find that Ur has an angle 
777,3 = 5r-tI\ 2 r a ^ er ^ ne fi rs t fold. Under condition (Q) we 
have mi > rri2 > and ms < 0. Therefore there must be 
a transverse homoclinic intersection in R. This implies an 
infinity of homoclinic intersections and the existence of a 
chaotic orbit. 

We now investigate the stability of this orbit. The basin 
boundary is formed by Sl- Sl folds at the y-axis and 
intersects the x-axis at point C. The portion of Ul to 
the left of L* goes to infinity and the portion to the right 
of L* leads to the chaotic orbit. Ul meets the x-axis at 
point D, and then undergoes repeated foldings leading to 
an intricately folded compact structure as shown in Fig.||. 

The unstable eigenvector at L* has a negative slope 
given by (— Sl/X\l). Therefore it must have a hetero- 
crine intersection with Sr. Since both Ul and Ur have 
transverse intersections with Sr, by the Lambda Lemma 
H we conclude that for each point q on Ur and for each 
e-neighborhood N e (q), there exist points of Ul in N e (q). 
Since Ul comes arbitrarily close to Ur, the attractor must 
span Ul in one side of the heteroclinic point. 

Since all initial conditions in L converge on Ul and all 
initial conditions in R converge on Ur, and since there 




Figure 3: The stable and unstable manifolds of L* for 
t l = 1.7, S L = 0.5, t r = -1.7, Sr = 0.5. R* is marked by 
the small cross inside the attractor. 

are points of Ul in every neighborhood of Ur, we con- 
clude that the attractor is unique. This chaotic attractor 
can not be destroyed by small changes in the parameters. 
Since small changes in the parameters can only cause small 
changes in the Lyapunov exponents, where the chaotic at- 
tractor is stable, it is also robust. 

It is clear from the above geometrical structure that no 
point of the attractor can be to the right of point D. If 
D lies towards the left of C, the chaotic orbit is stable. If 
D falls outside the basin of attraction, it is an unstable 
chaotic orbit or chaotic saddle. From this, the condition 
of stability of the chaotic attractor is obtained as 

^LTRXi L -SRXi L X2L J r^RX2L-^LTR J rT L S L -S 2 L -X2L^L > 

(5) 

If Sl = Sr = S this condition reduces to trX\l — X\lX2l + 

TL ~ Tr - S>0. 

The robust chaotic orbit continues to exist as tl is re- 
duced below (1+<5l)- With tl slightly below there 
is no fixed point in L for \i > but the invariant manifolds 
suffer only slight change. The invariant manifold of L as- 
sociated with Ail still forms the attractor. The invariant 
manifolds in L, however, cease to exist for tl < 2 ^/Sl since 
the eigenvalues become complex. As tl is reduced below 
2V^l there is a sudden reduction in the size of the attrac- 
tor as it spans only Ur. So long as Ul exists, multiple 
attractors can not exist and therefore if the main attractor 
is chaotic, it is also robust. 

Therefore we see that for 1>Jl>0, 1>Jr>0, the 
normal form (Q) exhibits robust chaos in a portion of pa- 
rameter space bounded by the conditions tr = — (1 + Sr), 
tl>2VSl and (|5|), as shown in Figj^. There is a symmet- 
ric region of the parameter space with the roles of R and 
L interchanged, where the same phenomena are observed 
for /i<0. 
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Figure 4: Schematic diagram of the parameter space re- 
gion of the normal form (||) where robust chaos is observed 
for 1> S L >0, 1> S R >0 and fi>0. 



For Sl < 0, the invariant manifolds Ul and Sl always exist 
as the eigenvalues are real for all tl- Therefore multiple 
attr actors can not exist for Sl<0. 

Since (|2|) is a normal form of the piecewise smooth map 
(0), it is expected that robust chaos would be observable 
in many piecewise smooth maps in the neighborhood of 
border collision bifurcations, provided that there are no 
more than one period- 1 fixed point in R a and R^, there 
exist homoclinic as well as heteroclinic intersections of the 
invariant manifolds associated with these fixed points, and 
the trace and determinant at the two sides of the border- 
line satisfy the above conditions. The example of the boost 
converter is a case in point. 

A major conclusion of this Letter is that one should 
use piecewise smooth systems in applications that require 
reliable operation under chaos. 



At very low values of the determinant, i.e., when the 
system is very close to being one-dimensional, the main 
attractor may not remain chaotic even for tl > 2\J~Sl, as 
periodic orbits become stable. The conditions for emer- 
gence of periodic windows for the one-dimensional case 
have been derived in ||, |To| . Therefore for one dimen- 
sional systems, the parameter range for robust chaos is 
bounded by tr = 1 , tr > — T T L L _ X , and the lower limit of tl 
is given by the conditions of existence of various periodic 
windows. Here the r's are to be interpreted as the slopes 
of the piecewise linear function in the two halves separated 
by x — 0. On the other hand, if the determinants in the 
two sides are unity, the region in the tl — tr space for 
robust chaos shrinks to zero area. 

The cases with negative determinant are investigated 
following the same method. For the sake of brevity, we 
present the results without much explanation. 

For — 1 < 5r < 0, we have 1 > Xir > 0, X2R < — 1, and 
R* is located above the x-axis. A positive value of Xir 
implies that XJl converges on XJr from one side. If 



tl-1-S l tr-1-Sr 

then the intersection of Ul with the x-axis remains the 
rightmost point of the attractor and (||) still gives the pa- 
rameter range for boundary crisis. But if (||) is not satis- 
fied, the intersection of XJr with the x-axis becomes the 
rightmost point of the attractor, and the condition of ex- 
istence of the chaotic attractor changes to 

X2R - 1 < $l (tl - S L - X 2 l) ^ 
tr-1-Sr (t l - 1 - Sl) (5rX 2 l - 5 l t r ) 

For Sl <0 and Sr < 0, L* is below the x-axis and the 
same logic as above applies. But if Sl < and Sr > 0, 
the stable manifold of R* has a negative eigenvalue and 
hence Ul does not approach XJr from one side. Therefore, 
if (H) is not satisfied, there is no analytic condition for 
boundary crisis — it has to be determined numerically. 
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